LIBRARY  USE  ONLY 


Ti-m-mi 

NUWC-NPT  TM  922089 

G>py  / 

NAVAL  UNDERSEA  WARFARE  CENTER  DIVISION 
NEWPORT,  RHODE  ISLAND 


Technical  Memorandum 


PROCEEDINGS  OF  THE 

NUWC  DIVISION  NEWPORT  SEMINAR  SERIES  ON 

TURBULENCE  AND  ITS  CONTROL 


1  October  1992 


Weapons  Technology  and 
Undersea  Systems  Dept. 


/ 


NAVAL  UNDERWATER  SYSTEMS  CENTER 
NEWPORT  LABORATORY 
NEWPORT.  RHODE  ISLAND  O284M047 
RETURN  TO:  TECHNICAL  LISRaSY 


Approved  for  public  release;  distribution  is  unlimited. 


LIBRARY  USE  ONLY 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

01  OCT  1992 


2.  REPORT  TYPE 

Technical  Memorandum 


3.  DATES  COVERED 

01-06-1992  to  31-08-1992 


4.  TITLE  AND  SUBTITLE 

Proceedings  of  the  NUWC  Division  Newport  Seminar  Series  on 
Turbulence  and  its  Control 

6.  AUTHOR(S) 

Promode  Bandyopadhyay;  J.  Meng;  A.  Hussain;  Steven  Orszag;  Daniel 
Nosenchuck 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Undersea  Warfare  Center  Division, Newport, RI, 02841 


5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

TM  922089 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSOR/MONITOR'S  ACRONYM(S) 

1 1.  SPONSOR/MONITOR'S  REPORT 
NUMBER!  S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

NUWC2015  Additional  authors:  Brown,  Garry;  Roshko,  Anatol 

14.  ABSTRACT 

This  memorandum  records  the  proceedings  of  a  four-part  seminar  series  on  turbulence  and  its  control 
sponsored  by  the  Naval  Undersea  Warfare  Center  Division,  Newport,  Rl,  during  the  summer  of  1992.  The 
Naval  Undersea  Warfare  Center  Division,  Newport,  RI,  organized  a  four-part  seminar  series  titled 
Turbulence  and  Its  Control  during  the  summer  of  1992  .  One  seminar  was  held  in  June,  one  in  July,  and 
two  in  August.  These  seminars  were  an  activity  of  NUWC’s  Hydrodynamics  Sphere  of  Excellence,  which  is 
one  of  the  Center’s  leadership  areas.  The  presentation  materials  used  during  the  seminars,  consisting 
mostly  of  informal  viewgraphs,  are  reproduced  in  this  report  in  their  original  form. 

15.  SUBJECT  TERMS 

Turbulence;  Hydrodynamics  Sphere  of  Excellence;  viscous  flow;  RNG;  renormalization  group  theory; 
magnetohydrodynamic  method;  bluff-body  wakes 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF 

OF  ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

205 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


ABSTRACT 
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FOREWORD 


The  Naval  Undersea  Warfare  Center  Division,  Newport,  RI ,  organized  a 
four-part  seminar  series  titled  Turbulence  and  Its  Control  during  the  summer  of 
1992.  One  seminar  was  held  in  June,  one  in  July,  and  two  in  August.  These 
seminars  were  an  activity  of  NUWC's  Hydrodynamics  Sphere  of  Excellence,  which  is 
one  of  the  Center's  leadership  areas.  The  presentation  materials  used  during 
the  seminars,  consisting  mostly  of  informal  viewgraphs,  are  reproduced  in  this 
report  in  their  original  form. 

In  the  first  seminar,  Professor  Hussain  discussed  in  great  depth  the 
viscous  flow  physics  of  the  numerically  simulated  "simple"  problem  of  the 
interaction  between  two  vortices. 

In  the  second  seminar,  Professor  Orzwag  presented  renormalization  group 
(RNG)  theory,  which  has  generated  a  great  deal  of  interest  recently  due  to  its 
surprising  success  in  the  rational  calculation  of  widely  different  turbulent 
flows.  RNG  theory  has  already  been  closely  scrutinized  by  many  independent 
groups.  In  high  Reynolds  number  practical  flows,  particularly  those  with 
strong  anisotropy,  RNG  might  become  the  standard  tool. 

The  third  seminar,  led  by  Professors  Nosenchuck  and  Brown,  focused  on  a 
novel  magnetohydrodynamic  method  of  turbulence  control.  Professor  Nosenchuck, 
with  his  keen  aptitude  for  application,  and  Professor  Brown,  with  his  deep 
insight,  made  a  formidable  team.  The  turbulence  control  they  have  achieved 
with  minimal  power  has  surpassed  that  obtained  with  polymers. 

In  the  fourth  and  final  seminar,  Professor  Roshko  reviewed  our 
understanding  of  the  turbulence  in  bluff-body  wakes  —  a  topic  that  has  always 
been  of  prime  interest  to  him.  He  discussed  his  recent  works,  emphasizing  the 
end-effects  in  laboratory  experiments  and  highlighting  the  lack  of  agreement 
between  so-called  two-dimensional  measurements  and  computer  simulations. 
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1.  New  Aspects  of  Vortex  Dynamics  and  Hydrodynamic  Turbulence 


A.  K.  M.  F.  Hussain 
University  of  Houston 
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NUWC  Division  Newport,  R.  I. 

SEMINAR  NOTICE 


NEW  ASPECTS  OF  VORTEX  DYNAMICS  AND 
HYDRODYNAMIC  TURBULENCE 


A.  K.  M.  F.  Hussain 
Cullen  Distinguished  Professor 
University  of  Houston 


We  try  to  shed  some  light  on  coherent  structures  and  turbulence  phenomena 
through  studies  of  the  new  aspects  of  vortex  dynamics,  and  of  coherent  vortex 
interaction  with  fine  scale  turbulence.  These  studies  are  done  by  direct  numerical 
simulation  of  the  Navier-Stokes  equations.  First,  we  explore  the  vortex  reconnec¬ 
tion  mechanism  and  its  role  in  turbulence  cascade  and  mixing.  We  come  to  realize 
that  core  dynamics  is  important  in  reconnection  and,  although  ignored  so  far,  is 
very  important  for  vortex  dynamics.  We  explain  core  dynamics  first  in  the  frame¬ 
work  of  traditional  quantities  as  colliding  wavepackets  resulting  from  coupling  of 
meridional  flow  and  swirl,  and  then  in  the  framework  of  a  new  mathematical  tool 
-  ’complex  helical  wave  decomposition’  -  which  gives  a  dearer  understanding  of 
the  flow  physics  in  terms  of  polarized  vortidty  waves,  expressed  in  terms  of  the 
eigenmodes  of  the  curl  operator.  Finally,  we  discuss  the  symbiotic  relationship 
between  coherent  structures  and  incoherent  turbulence  and  question  the  validity 
of  ’local  isotropy’  -  the  centerpiece  of  Kolmogorov’s  equilibrium  hypothesis  and 
of  virtually  all  theories  of  turbulence. 


Wednesday,  3  June  1992 
Conference  Room,  Bldg.  990  (6th  Floor) 

Time:  10:30  AM 

POC:  Dr.  Promode  R.  Bandyopadhyay  (Code  804;  x2588) 
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Figure  1,  Icol  surface  at  30%  of  the  initial  peak,  (a)  t=0; 
(b)  t=2.25;  (c)  t=3.5;  (d)  t=4.75;  (e)  t=6.0 


Figure  2.  Scalar  at  5%  of  the  initial  peak,  (a)  t=0; 

(b)  t=2.25;  (c)  t=3.5;  (d)  t=4.75;  (e)  t= 
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Figure  4.  kol  contours  in  Ka ;  contour  increments  Aco=4. 
(a)  t=3.0;  (b)  t=3.75;  (c)  t=4.5;  (d)  t=6.0 
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In  this  talk,  we  shall  describe  the  development  and  application  of  renormali¬ 
zation  group  (RNG)  turbulence  transport  models  with  non-equilibrium  rate-of- 
strain  effects  included.  A  wide  variety  of  examples  will  be  given,  including  tur¬ 
bulent  flows,  heat  transfer,  and  pre-mixed  combustion  in  step,  cylinder,  turn¬ 
around  bend,  valve,  constricted  pipe,  blunt  plate,  and  other  geometries.  The 
robustness  of  the  modeling  for  massively  separated,  non-equilibrium,  and  swirling 
flows  will  be  discussed.  The  relationship  of  transport  modeling  to  fUll  Navier- 
Stokes  solutions  and  large-eddy  simulations  will  also  be  discussed,  especially  in 
the  context  of  describing  coherent  structure  dynamics. 
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WORK  REQUIREMENTS  FOR  TURBULENCE  COMPUTATIONS 


Three-Dimensional  Turbulence  — 

Storage  0[(R3/4)3]  =  0[R9/4] 

Work  0[(R3/4)4]  =  0(R3] 

Two-Dimensional  "Turbulence"  — 

Storage  0[(R1/2)2]  =  0[R]  - 

Work  0[(R1/2)31  =  0[R3/2] 

GOAL  OF  TURBULENCE  THEORY: 

REDUCE  THESE  STORAGE  AND  WORK  REQUIREMENTS 


DIRECT  NUMERICAL  SIMULATION 

DNS 


A 

Grid 

size 


Navier- Stokes  Equations 


i  ~  It 
ld  ~  ju*FiL 
Dissipation 

scale 


Large  Scale 
(eddy) 
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CM2  Simulation  of  Homogeneous  Turbulence 

-Advanced  Computing  Laboratory- 
-Los  Alamos  National  Laboratory— 


Shiyi  Chen,  Xiaowen  Shan  LANL 

Robert  Kraichnan 

Zhen-Su  She,  Steven  Orszag  Princeton 


512  x  512  x  512  Spectral  Code  («  1+  min/timestep) 


Runs  to  date  (11/24/91):  R\  =  36,  70,  160 
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DRAG 


Figure  52:  Drag  on  each  waU  vs.  Reynolds  number  for  riblet  channel  simulation. 
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blet  valley 


SUBGRID  TURBULENCE 


Effects  of  Subgrid  Turbulence 

•Turbulence  Dissipation 
•Turbulence  Production  - 
Random  Forcing 


2-18 


LARGE-EDDY  SIMULATION 


LES 


Turbulence 
^ - Model  - 

4 


Solution  of 

LES-modified  Navier-Stokes 


L. 

Dissipation 

scale 


1 


A 

Grid 

Size 


J 


L 

Large 

Eddy 


Length-scale  model 
for  ‘subgrid ’  stresses 


supergrid  I 


Result:  U  pergnd  random  on  scales 

large  eddy 
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LARGE  EDDY  MODELS  OF  TURBULENCE 


—  —  Navier-Stokes  Equations  Solved  at  Large  (Super-Grid)  Scales 

—  One-Point  (Eddy  Viscosity,  K-£ Model)  Closure  for  Sub-Grid  Scales 

—  Only  Sub-Grid  Small  Scales  Removed  from  the  Dynamical  Equations 

—  Large  Numerical  Calculations  Involved  Solving  for  Super-Grid  Scales 

—  Applications  to  Homogeneous  and  Shear  Flows 

—  .Large  Scale  Structures  are  Computable  in  Detail 

—  Spectra  are  Computable  Up  to  Grid-Induced  Cutoff 
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TRANSPORT  MODEL 


REYNOLDS  AVERAGED  NAVIER-STOKES 


Turbulence 
<£ - Model  - 


l 


1 


I 


7  ~  1  T 

Dissipation 

scale 


« 


A 

Grid 

Size 


«  L 

Large  scale 
(eddy) 


Eddy  viscosity  or  K  —  S  model 
e.g.,  uv=  -ucddy Vt7 


£ - 

/ 


Large  eddies  solved 
by  Reynolds  averaged 
Navier-  Stokes,  *  . 


e-g 


dU_ 
»  8t 


dz  ^eddy) 


§U 

dz 


Result:  Smooth,  non-random  U 
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TRANSPORT  MODELS  OF  TURBULENCE 

—  Reynolds  Averaged  Equations  of  Motion 

—  Equations  for  Mean  Velocity,  RMS  Velocity  Fluctuations, .... 

—  Closures  Often  Based  on  Gradient  Transport  Ideas  (Eddy  Viscosity) 

—  All  Small  Scales  Removed  from  the  Dynamical  Equations 

—  Applications  to  Shear  Flows 

—  No  Information  Deduced  About  Small  Scale  Spectra, ... 
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ANALYTICAL  THEORIES  OF  TURBULENCE 


—  Multi-Point,  Muhi-Time  Moments  Usually  Involved 

—  All  Scales  Treated  Statistically 

i 

—  Renormalized  Perturbation  Methods  Normally  Used 

—  Currently  Applied  Principally  to  Homogeneous  Turbulence 

—  Huge  Numerical  Calculations  Involved  for  Shear  Flows 
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DYNAMIC  RENORMALIZATION  GROUP 


O  kQ  k  Ae'1  A 


Infrared  renormalization  group-  long-d;stance 
behavior 

1.  Remove  degrees  of  freedom 

Ae_i<q<  A 

2. Rewrite  (by  rescaling)  Navier-Stokes 

equations  as  renormalized  system  for 
v^^Ae*1)  with  modified  viscosity, 
force,  coupling 
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“The  renormalization  group  is  one  of  the 
fundamental  approaches  to  tackling  this 
problem  of  what  to  do  when  you  cannot 
make  your  grid  small  enough  to  use  the 
fundamental  equation.  How  do  you  in¬ 
crease  the  grid  spacing  beyond  the  level  of 
a  straight  numerical  approach ,  yet  preserve 
all  of  the  reliability  that  working  from  a 
fundamental  equation  can  give  you?” 

Kenneth  Wilson  (1985) 
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Intrinsic  Stirring  Force  in  Turbulence  and  the  e-Expansion 


^rg^scaje^cetomodenhe^ffect^f^/ffeTand^un^^conditionsJ 


(fa(k)f,(k'))«D0P^(k)5(k)5(k+k')5(ffl+ffl') 


DIFFICULTY:  Nonlinear  solutions  of  the  Navier-Stokes  equations  involve  an 

'infinity'  of  Interacting  f's  to  produce  u(k)  [k  finite] 

Correspondence  Principle  (Yakhot  &  Orszag  1 986)  | 


(#ll 

<5(k)  = 

■  lim (4  -PS''  1 

e-»4  4 7Z 

Gel'fand 


,  '  ..  iw  -*•'*. 


RENORMALIZATION  GROUP  (RNG)  THEORY 
-  A  Practical  Approach  to  the  Turbulence  Problem  - 
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Lange vin  Model 


dul 

dt 


-“v(e)”ue£~2  +  f 


Nonlocal  convolution  operator 


Assume 


dv(£) 

di 


oc  £ 


(rate  of  energy  dissipation) 


Dimensional  analysis 


Solution: 


dv  _  A££3 
~di  ~  I/2 


v(ld)  = 


1/3 


(i  >  id) 
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Renormalization  Group  Subgrid  Model  for  LES 

K  oc  1/A 


l*eddy  —  V mol 


,  3  AS 

1+  T 


4  (2t r) 


47/3 


A4- 


z/. 


Vd 


mol 


tti2 


1/3 


£  —  Z^ eddy  VU 


High  turbulence  limit 

Veddy 

Transitional  regimes 


3  A 


1  1/2 


4  (2tt)< 


A2  |VU| 


•&1,  =  J  (A?(A4  -  V >eddy  |VU|2 

— -  unphysical  roots  can  be  a  problem 


Alternate  Formulation 


S  = 


^eddy 


S  —  V eddy  VTJ 


4  _  3  , 

^ eddy  ^mol  ^ eddy  r  T 


3  AS2 


(A4  - 


4  (2ir)‘ 


>0 


Only  1  real  root  with  veddy  >  &w 
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2 


X 


Turbulent  flow  over  a  backward  facing  step 

Figure  3 


2-31 


APPLICATION  OF  RG  TO  TURBULENCE  MODELLING 


Prandtl  mixing  length  theory 


s  =  •'eddy  |VU|2 


k0  =  27r/L 


^eddy  — 


L  S.  ■  ^  / 


|vu| 


=  0.094 
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RNG-BASED  ALGEBRAIC  MODEL 


Eddy  Viscosity: 


//3 


au,  .  dUj  2  ,  „rr  eP 


e-"(a^  +  ax. 


(1) 

(2) 


Length  Scale  (A): 


where 


A  =  y+  if  y+  <  A,  and  =  A  otherwise , 


A = 7(1  -  b"'s' 

8.  -  displacement  thickness 
9  -  momentum-loss  thickness 
8./ 6  -  shape  factor 


(3) 

(4) 


Pressure  Gradient  Term  ( f3Uidp/dxi  ): 

/?  =  /?i(l.-exp(-£2(y+-A)/A))  (5) 


Constants: 


a  =  0.0256 
C  «  100 
7  =  0.28 
A  =5. 
A  =  0-12 
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.  friction  coefficient 


Experiments 


♦ 


_  Computation 

B>2dwln-Lomax  Model 


ARC  BAL  &  LOMAX 

MACH  0.675  ALPHA  0.000 
CL  -0.0701  CD  0.0156  CM  tUAH 
GRID  330X64  NCYC  300  W50.000E-00 


UtP  Of  O'  Co-  OSO  CO 


e  Experiments 

-  Computation 

PNG  A] (ebrelc  Model 


ARC 

MACH  0.675  ALPH  A  0.000 
CL  -0.0041  CD  0.0162  CM  0.1344 
GRID  330X64  NCYC  100  PX30.CWE.00 


Additional  Physics 


Example:  Stratified  Shear  Flow 


^ eddy 


=  VQ 


1/3 
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t  n-(5)n 


* 

3 


-  3 


-  4 


RNG  CLOSURE  FOR  STRATIFIED  TURBULENCE 


x  Theory  (<0;  •  Theory  C>0 

a  Experiment  ■  Experiment 


e  =  y/L 


I  Stable  Stratification 
I[  Unstable  Stratification 


.002  0.005  0.01 


2  ,  0.05  0. 

U 


APPLICATION  OF  RG  TO  TURBULENCE  MODELLING 


RG  Length  Scale  Model  for  High  Mach  Flows 

1C  =  k~2/3  ~  72c2  [V2  <  72c2  or  else  eddy  shoclets  form 

so 


so 


^eddy 


£1/3ko4/3~^ 


4  4 

*  <r 


T7I  2 


^eddy 


vu| 


^eddy 


rsj  • 


|vu| 


Self-Focusing  of  High  Ma  Jet /Wake  Flows 
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RNG  K-£  Eddy  Viscosity 
High  Turbulence  Formulation 


v 


(iAf)l7V3 


K  =  f°°  CKo^2/3k-5/3dk  =  -CKO£2/3i2/3 

Ji -i  2 


1/ 


c, 


(fA)1/3 

(§Cko)2 

(fA)1/3 

(|CKo)2 


K2 

1 

RNG 


0.0845 


General  Foundation 


v  —  C^K2 /£  only  if  v  i/q 


v  = 


i  +  */“ 


K 


"o  \/£ 


>  ^o) 
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/(enornicJiyJrtfn  (jrvup  (M&J  K'  £  Hohts 

((  -  -turbu/e^t  Icmefrt  |  Cr 

£  =  ■bM’bdt'd  ohrstptih**  <£=’>'  /W/ 

dO  , 

3)< 

-  ^  <  *  • 

CD-t 

2>£  __ 

^  2  •  •  % 

1>± 


%N(r  /fl>€ory  fix  £lx*Jitird  H»Jt  I  (€ty*ru.dl) 

{0nrA|M^r  wtWe«^ JL-^wipf 


-  ffcfi  o£j4r»  >h  -  i«*jp#f4a*W  «fif 

HOn*C|u«litrttf.iH  t^ftcti  Ck«(  rytrf  Arhfk 0*  Awi*  * 


-  ^3v*«4i*ftca.*fio^  ^  ro^«ihoir\  C'p'flcdz  4cco*u*UJ& 
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RNG  K -S  Model  with  Finite  Rate-of- Strain  Effects 


V.  Yakhot,  S.  Thangam,  T.  B.  Gatski, 

S.  A.  Orszag,  C.  G.  Speziale,  L.  M.  Smith 


dUi  „„  „  d 

—  4-  U  ■  V£7,  =  -Vip  + 

dK 


fdUj  dUj 

dxj  l"  \  dxj  "**  dxi 


di 


—  _  d  (  dK' 

+  U  •  VK  =  -S  -  TijSij  +  —  (  aKV- 


dxi 


dxr 


f  +  U  •  V?  =  -Cei  LrijSij -<?4-*+s;W 


Cei  =  1.42,  C£2  =  1.68,  C p  =  0.0845,  aK  =  cte  =  1.39 


7£  =  2i/oS{j 


dut  dui 
dxi  dxj 


Pade  Approximation  to  Expansion  of  TZ  in  Powers  of 


SK 


n=  ~  V/Vo) 

1  +  firf  K 
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Evaluation  of  9? 


•  Consistency  with  weakly  sheared  turbulence  77  — ►  0 

$ft  ~  i/S3 

•  Consistency  with  strongly  sheared  (rapid  distortion)  turbulence 

—  0(77)  77  — *  00 

Pade  Approximant 

R  =  ^tS3(1  -  77/770) 

1  +  Pv2 

(Y,rnr]n  is  a  geometrical  series) 

»=  C/xT73(1  —  77/770)  ^ 

1  +  /3t73  K 
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RNG  FORMULATION  OF  K-e  TURBULENCE  MODEL 


de/d t  +  Uk3£/9xk  =  a3Uj/9xk  <u;uk>  -Y  +  9/9xk  a£v  9e/9xk 

9k/9t  +  Uk9k/9xk  =  9Uj/9xk  <u;uk>  -  £  +  9/9xk  oc^v  9k/9xk 

d(a He)  =  -0.2176  dv/Y-j 
d(Y/Ve3)  =  -0.3089  dv/Y1 
d(  k/Ve)  =  1.63vdv/Y-| 

.Y^Kv/v^-I+C]1^ 

|(ak-  1.3930)/0.3930|0-6321 
•  |(ak+  2.3930)/3.3930|0-3679  =  vQ/vR 
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K*(t*) 


5 

4 

3 

2 

1 

0 

0  2  4  6  8  10 

t* 

Time  evolution  of  the  turbulent  kinetic  energy  In  homogeneous  shear  flow. 
- Relaxation  model;  o  Large-eddy  simulation  of  Bardina  etalu  for  e0/S/C0= 0.296 

Figure  2 
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BACKWARD-FACING  STEP: 


Y/H 


-2  0 


5 


1  0 


X/H 

(a)  Streamlines 


X/H  :  1.33  2.67  5.33  6.22  8.00  10.67 


01  01  01  Ol  01  01 


j; 

Uo 

(b)  Dimensionless  mean  velocity  profile 

( - Computations  with  Isotropic  eddy  viscosity; 

o  Experiments  of  Kim  ef  al,  1980;  Eaton  &  Johnston,  1981) 


Computed  mean  flowfleld  for  the  new  RNG  K-z  model 
[E  =  1:3;  Re  =  132,000;  200x100  mesh] 

Figure  4 
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FlGL'KE  10.  Aluminum-powder  pictures  of  streamlines  over  a  step.  Exposure  time  0.5  (upper)  and  5  seconds  (lower). 


(a)  Reattachment  at  XH  -  “  J 


(b)  Primary  recirculation  zone 


(a)  Corner  eddy 

Streamlines  for  the  flow  over  a  backward  facing  step  computed 
using  the  RNG  K-  8  model  (Re  =  88000,  E=2:3) 
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BACKWARD-FACING  STEP: 


X/H  :  1.33  5.33  7.57  8.55  10.33  15.66 


(a)  Turbulence  Intensity 


X/H  :  1.33  5.33  7.67  8.55  10.33  15.66 


0.0  0.01  0.02  0.01  0.02  0.01  0.02  0.01  0.02  0.01  0.02  0.01  0.02 


(b)  Turbulence  shear  stress 

Computed  turbulence  slressesforthenewRNGK-F.modeI[E= 1:3;  Re  =  132,000; 

200x100  mesh;  -  computations  with  anisotropic  eddy  viscosity; 

o  experiments  of  Kim  et  al,  1980;  Eaton  &  Johnston,  1981] 

Figure  7 
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h*/hfi,hf/h 


RNG  K— e  model  with  RNG  heat  transfer  for  air,  cr0  =  0.7) 


Heatflux  ratios  for  the  air  flow  over  a  backward  facing  step  computed  using 
the  RNG  K-e  model,  Re  =  88,000,  E  =  2:3,  a  =  0.7.  (hj3  =  heat  flux 
at  the  step  corner,  hje  =  heat  flux  at  the  exit  ) 
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hf/h{,,h(/h{,  h(/h,,,hf/h 


RNG  K — c  model  with  RNG  heat  transfer  (liquid  metal,  a0  =  0.025) 


Standard  K— e  model  with  standard  heat  transfer 
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'  '  I  ■  '  , '  :  .  ■  ■ 

v  ‘  v  V,  ■  ‘ 

■ 7  ■  iO"'  ■  ■  v 

Streamlines  for  the  flow  in  a  constricted  pipe  computed  using  the 
standard  K~  F  model  (  Re  *  10,000,  I)  throat  / 1)  jn  =  0.5) 


Streamlines  for  the  flow  in  a  constricted  pipe  at  computed  using  the 
RNC.K-f  Model  (Re -10,000 , ,  =  0.5) 
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Streamlines  for  the  flow  in  a  turnaround  duct  computed 
i  using  the  standard  K- E  model 
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Streamlines  tor  the  flow  oven  a  blunt  flat  plate  computed*  using 
a)  theRNGK  .£  model  b)  the  standard  K  £  model  (Rc  -  21,600) 


Streamlines  in  the  wake  of  a  cylinder  at  Re~l  4,500  computed 
_ _  using  the  RNG  K-  £  Vtodei 
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Intensity  (  <uu>1/2/U( 


X/D 

Comparison  between  the  RNG  and  standard  K-e  models  of  the  average 
turbulence  intensity  in  the  wake  of  the  cylinder  at  Re  =  14,500. 
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•  .  ,  .  ,v 

Pressure  contours  in  the  flow  over  a ’cylinder  at  Re  =  i4,5U0  computed 

using  the  RNG  K  el  Model 


Fig.  ii  The  computed  and  experimental  flame-  Fig.  12  The  computed  and  experimental  flame- 
propagation  process  in  a  4-valve  engine  propagation  process  in  a  4-valve  engine 

(CASE  A-ll)  Engine  speed  =1400rpm,  (CASE  A-lll)  Engine  speed  =21  OOrpm 
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Crank  angle  al  the  maximum  pressure 

Fig. 14  The  comparison  of  the  computed  maximum 
pressure  with  the  experimental  data 
(Engine  speed  variation) 


700 rpm 
(CASE  A-l) 


□  Computation 
I  Experiment 


1400rpm 
(CASE  A  ll) 


2100rpm 
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Crank  angle  at  the  maximum  pressure 


Fig.  16  The  comparison  of  the  computed  maximum 
pressure  with  the  experimental  data 
(A/F  ratio  variation) 


Fig.  15  The  comparison  of  the  computed  maximum 
pressure  with  the  experimental  data 
(Ignition  timing  variation) 
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Ignition  timing  variation 


340  360  300  400  420  440 

Crank  angle  degATDC 


Enoinen  speed  variation 


A/F  ratio  variation 


“■  Computed  pressure  (1st  cycle) 

—  Computed  pressure  (2nd  cycle) 

—  Experiments  pressure 

■«.  Computed  heat  release  rate  (1st  cycle) 
~  Computed  heat  release  rate  (2nd  cycle) 
Experimental  heat  release  rate 


Fig.  13  The  time  histories  of  the  computational-  pressure  and  heat  release  rate 
and  the  corresponding  experiments 
(a)  CASE  A-l.  (b)  CASE  A-ll,  (c)  CASE  A-lll, 

(d)  CASE  B,  (e)  CASE  C-l,  (f)  CASE  C-ll 
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Reverse  Blank 


3.  Direct  Control  of  Wall  Shear-Stress  in  a  Turbulent  Boundary  Layer 

Daniel  Nosenchuck  and  Garry  Brown 
Princeton  University 
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NUWC  Division  Newport,  R.  1. 

SEMINAR  NOTICE 

DIRECT  CONTROL  OF  WALL  TURBULENCE  USING 
ELECTROMAGNETIC  FORCES 


Professor  Daniel  M.  Nosenchuck 
and 

Professor  Garry  L.  Brown 
Princeton  University 

A  new  concept  and  technique  has  been  developed  to  directly  control  turbu¬ 
lence  production  in  boundary  layers.  Near-wall  vertical  motions  are  directly 
suppressed  through  the  application  of  a  Lorentz  force.  Current  (j)  and  mag¬ 
netic  ( B )  fields  parallel  to  the  boundary  and  normal  to  each  other  produce 
a  finite  Lorentz  force  j  X  B  normal  to  the  boundary.  Experiments  have 
been  performed  on  flat-plate  turbulent  boundary  layers  at  Ree  =  1700. 

With  the  application  of  modest  field  densities  (eg.  \B\  <  500  gauss  and 
\j\  <  10  mA/cm2),  measured  reductions  in  turbulent  stresses  within  the  con¬ 
trol  region  are  seen  to  exceed  90%.  Laser-sheet  flow  visualization  confirms 
the  substantial  reductions  in  turbulent  motion  at  y+  £  15.  It  is  suggested 
that  the  principal  reason  for  the  observed  effects  is  due  to  the  direct  control  of 
the  coherent  motions  responsible  for  turbulence  production  in  the  near-wall 
region. 


Thursday,  24th  September  1992 
Conference  Room  A,  Bldg .  102 
Time;  10:30  AM 


POC:  Dr.  Promode  R.  Bandyopadhyay  (Code  8234;  x2588) 
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Overview 


Turbulence  is  characterized  by: 

—  periodic  eruptions  of  unstable,  low-momentum 
‘near-wall’  fluid 

—  subsequent  inrush  of  high-speed  ‘outer-flow’  fluid 

—  resultant  large  skin-friction  drag 

Lorentz  force  easily  generated: 

—  surface  electrodes  produce  electric  field  with 
current  density  j 

—  magnetic  field  B  is  generated  parallel  to  surface 
and  normal  to  electric  field 

—  resultant  normal  force  is  j  x  B 

Direct  application  of  wall-normal  force  could 
prohibit  lift-up  and  bursting  of  near-wall  fluid 


Li**  Jl 

S/tcbro  \fh\ 
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McDonnell  Douglas 
Technologies  Incorporated 


TFM  Detailed  Concept 


Princeton 

University 


Elementary  Concepts 

MHD  Similarity 


Assume  1-D  flow;  applied  magnetic  field  Bo,  char¬ 
acteristic  length  L,  velocity  U ,  conductivity  a,  den¬ 
sity  p,  viscosity  p,  permeability  C ^ 


Per  unit  volume: 

Fv (viscous  force) 

^(inertia  force) 


U 

L2 

U2 


Fem  (electromagnetic  force)  <T  BqU 
Magnetic  Reynolds  Number  ( Rem ) 

induced  magnetic  field)  __  C^aBoUL 


— 


applied  magnetic  field 


Bo 

=  C^aUL 


Hartmann  number  (Ha) 


Ha  = 


em 


i  fv 


=  B0L 


a 


W 


Interaction  Parameter  (!)  (or  Stuart  Number) 


cB\L 

PU 


Theoretical  Considerations 


•  MOMENTUM 


— Vj>  +  L  +  /zV2u 


where  L  =  j  x  B  (Lorentz  force) 


•  VORTICITY 

Duj 


P~jy^  —  Pu  *  Vu  +  V  x  t  +  /iV  w 


If  L  =  Z^e^: 


.DcJt  .  .  ,  o 

-  =  pa;'-  Vu  -  — 1  +  pV2wx 


Pi 


dz 
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Turbulent  Boundary  Layer 
Energy  Balance 


////////////////////  :  MTC  Control  Layer 
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Theoretical  Considerations  - 

Continued 


ENERGY 

Near  the  wall 

,  ,du 

—u  v  — 

8y 


L'yv' 


+ 


Production  ~  Lorentz  Work  +  Dissipation 


this  suggests  the  nondimensionaf  parameter 

A  L  v 

P  ur 

for  scaling  Magnetic  Turbulence  Control 
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Some  Engineering  Considerations 

•  Fluid  must  be  electrically  conducting: 

—  sea-water  is  an  ideal  electrolyte 

—  small  amounts  of  electrolytes  injected  near  wall 
in  fresh-water  experiments 

—  gaseous  flows  may  be  seeded  with  ions  if  plasma 
is  not  naturally  present 

•  MTC  power  infinitesimal  since 

—  little  mechanical  work  is  done:  j  x  B  •  u  ~  0 

—  Joule  heating  (ere2)  can  be  made  negligibly  small 

•  j  and  B  fields  may  be  DC  or  AC 

•  Distributed  electrodes  and  magnet  poles  could  gen¬ 
erate  fields  that  conform  to  non-planar  surfaces 
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Experimental  Investigation  of  MTC 


•  Turbulent  boundary  layer  and  turbulent-spot 
flows  were  studied  on  a  flat  pJate  in  water 

•  Control  region  8  cm  wide  and  16  cm  long 

500  x  1,000  wall  units)  on  (£  of  plate 
200  cm  from  L.E. 

•  Control  zone  consisted  of  a  lucite  plate  with 

—  permanent  magnets  =>  transverse  field 

\By  ~  500  gauss) 

—  stainless-steel  surface-mounted  electrodes 

=>  longitudinal  field  (0.003  <  jx  <  20  mA/cm2) 

•  800  <  Ree  <  1700  upstream  of  contrci  zone 

•  Dilute  HCI,  NaCI,  and  NaOH/fluorescein  electrolyte 
injected  normal  to  wall  upstream  of  control  zone 

•  Diagnostics  included  hot-film  probes,  gauss-meters, 
electric  field  sensors,  laser-sheet  flow  visualization 
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MTC  Zone  on  a  Flat  Plate 


Electrolyte  Anode  Control  Region  Cathode 


a:  Planform  View 

(«  =  Jv  B) 
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Key  Results 


30  40  50 

Time  in  Seconds 


60  70  80 


a:  Control  Sequence:  Off-On-Off 


0  5  10  15  20  25  30  35 

Time  in  Seconds 


b:  Control  Sequence:  Off- On- Reverse 

MTC  on  Centerline  at  t/+  ~  1 
Velocity  vs  Time 
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x:  No  Electrolyte 
+:  Electrolyte  On;  No  Control 
□  :  Control  On:  j  x  B  <  0 


a:  UT  Computed  from  Baseline  Profile 


b:  uT  Computed  for  Each  Case  Separately 

Effect  of  MTC  on  Mean  Turbulent  Boundary  Layer  Velocity 

Profile 
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0  ‘ - - - i - ' - * - 1 - _ 

200  400  600  800  1000  1200 

t+ 


Effect  of  Magnitude  of  Lorentz  Force 
on  Magnetic  Turbulence  Control 

Velocity  vs  Time 

y+  ~  |  Reg  ~  1700 
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0  10  20  30  40  50 

Time  in  Seconds 

Turbulent  Motions  and  Stresses 

Velocity  vs  Time 
y+  ~  3  Ree  a  1100 


100 


Direct  Magnetic  Control 

of  a  Turbulent  Boundary  Layer 

^  y  ^  Planform  View  ^  y  ^  ^ 

y+<  5 


Elevation  View 

+ 

y 

200 


Laser-sheet  Flow  Visualization 
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Control  Zone 


o 

U 

^r(ref) 

0 

0 

o 


-40  -20  0  20  40  60  80  100 

x  (cm) 

Distance  from  trailing  edge  of  rear  electrode 


MTC  Relaxation 

Near-wall  Longitudinal  Velocity  (Wall-Shear)  Distribution 

On  Control  Zone  Centerline 

(Preliminary  Results) 


3-20 


Other  Applications 


Separation  control 

—  decrease  or  increase  form  drag 

—  provide  directional  control  through  net  lateral 
forces  and  moments 

Acoustic  field  attenuation/enhancement/tailoring 
Modification  of  unsteady  surface  pressure  fields 
Turbulent  fluid  mixing 
Surface  heat-transfer  control: 

—  j  x  B  <  0  =>  decreased  wall  heating 

—  j  x  B  >  0  =>  enhanced  wall  heating 


Summary  of  MTC 


•  A  new  concept  and  technique  using  electro-magnetic 
fields  has  been  developed  to  suppress  wall-turbulence 


•  Non-intrusive  magnets  and  surface  electrodes: 

—  impart  normal  force  j  x  B  near  wall  in  turbulent 
boundary  layers 

—  turbulence  generation  directly  suppressed 


•  Proof-of-concept  experiments  have  been  performed 
in  a  water  channel  on  a  flat  plate  model 

—  >  90%  reduction  in  turbulent  skin  friction  (drag) 

—  near  wall  velocity  fluctuations  attenuated  >  5x 

—  flow  visualization  indicates  nearly  ‘laminar’  flow 

—  minimalpower/current  required 
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Questions 


•  High  Reynolds  number  behavior - 

•  Dependence  on  ion  concentration  and  uniformity 


Underlying  dimensionless  parameters 


•  Relaxation ,  tiling,  and  three-dimensionality 

•  •  , 

_ _ _ _ _ ,  i 

•  Frequency  response; -field  attenuation - 

— •  Numerical  simulation - - - 
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McDonnell  Douglas 

Technologies  Incorporated 


TFM  Application  Concept  -  Sea 


Princeton  * 
University 


•  SAIL  FLOW 
SMOOTHING 


•  PROP/CONI  ROL 
FLOW 


SONAR  FLOW 
QUIETING 


4.  Bluff-Body  and  Other  Shear  Flows 


Anatol  Roshko 

California  Institute  of  Technology 
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SEMINAR  NOTICE 


BLUFF-BODY  AND  OTHER  SHEAR  FLOWS 


Anatol  Roshko 

Theodore  von  Kdrmdn  Professor  of  Aeronautics 

California  Institute  of  Technology 


The  old  problems  of  the  far  wake,  near  wake  and  forces  on  bluff  bodies  will 
be  discussed,  mainly  the  so-called  two-dimensional  flows  behind  circular  cylinders 
and  flat  plates.  Progress  and  problems  in  understanding  the  dynamics  of  the  flow 
will  be  discussed.  Topics  include:  effects  of  Reynolds  number;  effects  of  three 
dimensionality  in  nominally  ’two  dimensional  flows’;  roles  of  vortex  dynamics  and 
Reynolds  stress;  insights  from  experimental-computational  interactions  and  from 
other  shear  flows. 


Monday,  28th  September  1992 
Conference  Room  B,  Bldg.  102 
Time:  10:30  AM 


POC:  Dr.  Promode  R.  Bandyopadhyay  (Code  8234;  x2588) 
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FLOW  ABOUT  A  2-D  FLAT  PLATE 


-Carriage 


Water  surface7 


-  Laser  sheet 


U  =  aTa 


Flat  plate 


Tank  floor 


Side  view 


U - -j 


T*  =T 

a  -*•  a. 


Reu  =  — 

V 


Re  n  = 


l\fai 


tn  =  t- 


Tbp  yj 


view 


Motion  time  history 


Laser  sheet 


Rea  =  3966 
Rey  =  6096 
T*  =  1.54 
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END  EFFECTS  ON  THE  FLOW  PAST  FLAT  PLATES 


Water  surface 


Water  surface 


Water  surface 
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Comparison  of  Experiment  and  Numerical  Model 
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-o.S  ©  o.C 


-0.5 


0.5 


Pl*£ 


Force  Histories  On  A  Normal  Blunt  Plate 


After  Start-Up  In  Towing  Tank;  Re  =  3600  "Steady-State"  In  Free  Surface  Water  Tunnel;  Re  =  3060 


Position  x/C  x/C.  =  Ut/C 


T.  =  Strouhal  Period 


Force  Fluctuation  -  Free  Surface  Water  Tunnel 
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Histogram  of  Drag  Coefficient  Fluctuations  in  FSWT  -  AoA  =  90 


1.6 


uia/# 


lations  -  AoA  =  87.5  deg 


Numerical -Model- 
18000  Points 
Cd  =  3.19 


Density  Stiratlfijcation  Experiments 


•  reduce  3-d  turbulence  by  damping  span  wise 
motions. 


dp 

g—  =  constant  (spanwise  homogeneous) 

(Jj/C 


Richardson  Number:  measure  of  effect  of  bttayancy 

R  =  9  ^  =  n  Ap  gd  d 

1  P&  Pf  P  u2h 

d  _  A  pg  <P 
1  p  hU2 


^  (- 


4-21 


/ 


4-22 


Stratification  Effects  (Small  End  Gap, 


o 


20 


Position  X/C 


=  5cm,  AoA  =  90.0  deg) 


=  0.00,  Re  =  1000 
=  0.00,  Re  =  2500 
=  0.80,  Re  =  2500 
=  1.35,  Re  =  2500 
=  5.01,  Re  =  1000 
=8.45,  Re  =  1000 


<t  —  .5  Re  =  3300 


a  =  .5  Re  =  3300 
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IMPLICATIONS  AND  USES  OF  COHERENT  STRUCTURES 
(ideas  not  accessible  from  traditional  approaches) 

•  Role  of  global  (primary)  instability 

characteristic  for  each  shear  flow 
(implications  for  universal  Re  models?) 

•  Rational  explanation  and  correlation  of  parametric  effects: 

velocity  ratio  U2/Ui 
density  ratio  p^fpi 
compressibility  Mi ,  M2 
entrainment  parameters 

•  Sensitivity  to  perturbation  (¥fygnanski  and  Oster) 

uniqueness 

self-excitation?  (Corcos  and  Kaul) 

•  Controllability 

•  Lagrangian  view  useful  or  necessary 

e.g.  coherent  structure  as  a  stirred  reactor 

•  s 

(Broadwell) 

•  Closed  models  (Morris,  Liu) 

with  pu'v'  from  instability  equations 
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55  Invited  lectures  at  international  meetings  (including  keynote  lectures) 

Invited  seminars  (over  150  in  USA  and  abroad) 

59  Competitive  research  grants 
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STEVEN  A.  ORSZAG 


Steven  A.  Orszag  is  the  Hamrick  Professor  of  Engineering  and  Director  and  Professor  of 
Applied  and  Computational  Mathematics  at  Princeton  University.  He  studied  at  M.I.T.  and 
Cambridge  University  prior  to  receiving  his  Ph.D.  in  Astrophysics  from  Princeton  in  1966. 
After  a  year  at  the  Institute  for  Advanced  Study,  he  returned  to  M.I.T.  where  he  was 
Professor  of  Applied  Mathematics  until  assuming  his  present  position  in  1984.  His 
research  interests  include  numerical  analysis,  applied  mathematics,  and  fluid  dynamics.  His 
major  contributions  include  the  development  of  spectral  numerical  methods,  the  theoretical 
analysis  of  the  mechanism  of  transition  in  shear  flows,  the  first  numerical  simulations  of 
three-dimensional  turbulent  flows,  and  the  development  of  dynamic  renormalization  group 
methods  for  turbulence.  His  recent  awards  include  the  1986  AIAA  Fluids  and 
Plasmadynamics  Award,  a  Guggenheim  Fellowship,  and  the  1991  Otto  Laporte  Award  of 
the  American  Physical  Society.  Prof.  Orszag  is  Editor-in-Chief  of  the  Journal  of  Scientific 
Computing,  Amer.  Inst,  of  Physics  Series  in  Computational  Physics,  and  Springer  Series 
in  Computational  Physics.  He  has  written  over  250  papers  and  9  books.  He  is  co-author 
with  Carl  M.  Bender  of  the  widely  used  Advanced  Mathematical  Methods  for  Scientists 
and  Engineers  and  the  forthcoming  Partial  Differential  Equations  for  Scientists  and 
Engineers. 
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DANIEL  M.  NOSENCHUCK 


(office):  (609)258-5136 


[Pll  Redacted] 


EDUCATION 

Ph.D.  California  Institute  of  Technology,  Pasadena,  California 

Aeronautics,  Thesis  Advisor:  H.W.  Liepmann  /  June  1982 
M.S.  California  Institute  of  Technology,  Aeronautics  /  June  1977 
B.S.  Syracuse  University,  Syracuse  New  York 

Aeronautical  and  Mechanical  Engineering,  cum  laude  /  May  1976 

EXPERIENCE 

7/88  -  Present  Associate  Professor  of  Mechanical  and  Aerospace  Engineering 
8/83  -  6/88  Assistant  Professor 

Princeton  University,  Princeton  New  Jersey 

Teach  graduate  and  undergraduate  courses  in  fluid  mechanics. 
Current  research  areas: 

1)  development  and  application  of  the  Navier-Stokes  Computer 

2)  experimental  and  numerical  fluid  mechanics: 

control  of  complex  turbulent  and  vortical  flows 


HONORS 

-  Presidential  Young  Investigator  Award  (NSF),  1984-1989 

-  GTE  Emerging  Scholar,  1987 

-  Rheinstine  Award,  School  of  Engineering,  Princeton  University,  1986 

-  IBM  Faculty  Development  Award  1984  -  1985 

-  EMMY  Award  from  the  Academy  of  Television  Arts  and  Sciences  for 
Outstanding  Individual  Achievement  -  Special  Visual  Effects: 

‘ The  Day  After,'  1984 

-  The  William  F.  Ballhaus  Prize, 

for  an  ‘Outstanding  Doctoral  Dissertation  in  Aeronautics,’ 

California  Institute  of  Technology,  1982 

-  Graduate  Fellowships,  California  Institute  of  Technology,  1976 
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MISCELLANEOUS 


Patents: 

U.S.  Patent  No.  4811214  (March  1989):  ‘Dynamic  Reconfiguration  System  for 
Pipelined  Computers’,  Numerous  foreign  patents  also  awarded. 

U.S.  Patent  Pending:  ‘A  Parameterized  Optimizing  Compiler’ 

Reviewer  for: 

-  AIAA  Journal 

-  Journal  of  Fluid  Mechanics 

-  Physics  of  Fluids 

-  Review  of  Scientific  Instruments 

Department  of  Defense  Activities: 

-  Member  of  Defense  Science  Board  Ballistic  Missile  Defense  Task  Force 

-  Member  of  Defense  Science  Study  Group 

Consultant  to: 

-  Defense  Science  Board,  Washington  DC 

-  DSV  Partners,  Princeton,  NJ 

-  Institute  for  Defense  Analysis,  Alexandria  VA 

-  Northwest  Research  Associates  Inc.,  Belleview  WA 

-  Praxis  Film  Works,  North  Hollywood,  CA 

-  Union  Camp  Corporation,  Lawrenceville,  NJ 

Transferred  Navier-Stokes  Computer  technology  from  Princeton  University  to  Con¬ 
current  Computer  Corporation,  Tinton  Falls,  NJ,  in  1987,  and  Supercomputer 
Solutions  Inc,  San  Diego,  CA.,  in  1989. 
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DANIEL  M.  NOSENCHUCK 


Associate  Professor 

Daniel  M.  Nosenchuck’s  research  interests  include  experimental  and  computational  fluid 
mechanics,  dynamic  flow  visualization,  and  advanced  supercomputer  architectures.  His  major 
contributions  include  experimental  active  laminar-flow  control,  the  first  successful  demonstra¬ 
tion  of  active  turbulence  control,  and  the  development  of  a  parallel-processing  supercomputer. 
He  was  a  charter-year  recipient  of  the  five-year  NSF  Presidential  young  Investigator  Award 
(1984-89).  He  also  received  of  the  IBM  Faculty  De\-elopment  Award  (1984-85)  and  the  Prince¬ 
ton  University  School  of  Engineering  Rhinestein  Award  (1986)  for  work  related  to  the  im¬ 
plementation  of  unique  flow  fields,  and  the  development  of  new  flow  visualization  techniques. 
He  received  the  National  EMMY  Award  in  1984  for  Outstanding  Individual  Achievement  in 
Special  Visual  Effects.  He  is  active  in  industrial  consulting  and  also  consults  directly  with  the 
Department  of  Defense,  and  is  a  member  of  DoD  Defense  Science  Board  Task  Forces. 

The  underlying  theme  of  his  work  revolves  around  the  control  of  complex  fluid  flows.  To 
achieve  this,  he  is  currently  engaged  in  several  areas  of  research.  These  include  the  study 
and  experimental  active  control  of  turbulent  boundary-layers  in  a  low-speed  water  channel, 
wake-vortex  prediction  and  control,  the  development  of  a  new  three-dimensional  dynamic 
flow- visualization  technique,  and  the  design  and  construction  of  a  very-high-speed  computer 
for  use  in  complex  flow  simulations  and  control  applications. 

Experimental  Fluid  Mechanics 

The  goal  of  this  research  is  to  reduce  skin-friction  drag  by  modifying  turbulent  boundary 
layers.  A  new  technique  using  a  wall-normal  electromagnetic  body  force  to  directly  suppress 
near- wall  turbulent  instabilities  is  in  the  early  stages  of  development.  Preliminary  experimental 
.  results  indicate  the  possibility  of  dramatic  reductions  in  turbulence  and  wall-shear  with  little 
power  expenditure.  Other  experiments  involve  arrays  of  thin-film,  sensors  and  actuators  cou¬ 
pled  to  high-speed  feedback-control  electronics.  With  appropriate  special-purpose  electronics, 
a  significant  attenuation  of  boundary-layer  instabilities  has  been  realized.  A  second  major  area 
of  investigation  is  the  mitigation  of  the  wake- vortex  hazard  produced  behind  large  aircraft,  and 
the  alleviation  of  trailing  vortices  behind  submarines.  This  is  being  investigated  using  a  new 
approach  to  induce  instabilities  in  the  vortex  sheet  produced  by  lifting  surfaces.  Through  the 
development  and  application  of  an  optimizing  numerical  design  and  flow  simulation  procedure, 
small  modifications  to  the  lifting  surface  geometry  were  predicted  to  modify  substantially  the 
trailing  vortex.  Preliminary  experiments  were  used  to  validate  the  predictions. 

To  understand  these  complex,  nonsteady  boundary-layer  and  vortex  flows,  a  new  diagnostic 
technique  was  developed  within  the  lab.  A  laser-sheet  is  rapidly  scanned  through  flow-fields 
into  which  a  laser-fluorescing  dye  has  been  added,  and  a  series  of  two-dimensional  images  are 
obtained.  A  single  scan  is  comprised  of  many  closely-spaced  sheets.  From  this  space-filling 
data  set,  a  three-dimensional  flow  visualization  image,  representative  of  one  instant  in  the  flow, 
is  obtained.  Repeated  scanning  is  used  to  create  a  complete  three-dimensional  nonsteady  qual¬ 
itative  record  of  the  flow.  He  is  currently  extending  the  technique  to  encompass  quantitative 
methods,  and  investigating  its  use  as  a  real-time  input  sensor  for  active  control  experiments. 
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Navier-Stokes  Computer 

Problems  in  fluid  mechanics  involving  complex  flow  simulations  require  far  more  speed  and 
capacity  than  that  provided  by  current  and  proposed  conventional  and  parallel  supercomput¬ 
ers.  To  address  this  concern,  the  Navier-Stokes  Computer  (NSC)  was  developed.  The  NSC 
is  a  fully  general-purpose  parallel-processing  machine,  comprised  of  individual  Nodes,  each 
comparable  in  performance  to  current  supercomputers.  The  projected  speed  and  capacity  of 
a  128-Node  NSC  is  many  orders-of-magnitude  greater  than  that  of  existing  supercomputers. 
New  architectural  features  have  provided  the  capability  to  efficiently  address  a  far  greater 
range  of  problems  than  possible  on  conventional  machines.  The  first  large-scale  applications 
of  the  NSC  involved  the  simulation  of  complex  flows  at  moderate  Reynolds-numbers,  includ¬ 
ing  nonsteady  flows  with  separation  over  a  large  domain.  A  fast  optimizing  NSC  FORTRAN 
compiler,  which  uses  a  new  ‘approximate  simulation’  technique  coupled  is  in  development. 

Selected  Publications 

“The  Direct  Control  of  Wall  Shear  Stress  in  a  Turbulent  Boundary  Layer,”  (with  G.L. Brown), 
submitted  to  Phys.  Rev.  Lett.,  1992. 

“Parameterized  Memory/Processor  Optimizing  FORTRAN  Compiler  for  Parallel  Computers,” 

To  be  published  in  Proceedings  of  1992  Scalable  High  Performance  Computing  Conference, 
IEEE  Computer  Society  Press. 

“Control  of  Wingtip  Vortices,”  (with  W.S.  Flannery  and  G.L.  Brown)  Proceedings  of  FAA 
Wake  Vortex  Symposium,  Washington,  DC,  29-31  October  1991. 

“The  Navier-Stokes  Computer,”  Special-Purpose  Computers,  (with  W.S.  Flannery  and  E.  Hay- 
der)  B.  Alder,  ed.,  pp  97  -  134,  1988. 

“Active  Control  of  Sublayer  Disturbances  Using  an  Array  of  Heating  Elements,”  (with  M.  K.  Lynch 
and  J.  P.  Stratton)  Proc.  2nd  ASME/JSME  Thermal  Engineering  Joint  Conf.,  1987. 

“Three-Dimensional  Flow  Visualization  Using  Laser-Sheet  Scanning,”  (with  M.  K.  Lynch). 
Proceedings  of  the  AGARD  Conference  on  Aerodynamic  and  Related  Hydrodynamic  Studies 
Using  Water  Facilities,  AGARD  Conference  Preprint  No.  413,  pp  18-1  -  13,  1986. 

“Active  Control  of  Laminar-Turbulent  Transition,”  (with  H.W.  Liepmann).  J.  Fluid  Mech. 
Vol  118,  pp  201-204,  1982. 


A-10 


GARRY  L.  BROWN 

Professor  and  Chairman 


Garry  L.  Brown  was  appointed  Chairman  of  the  Department  of  Mechanical  & 
Aerospace  Engineering  at  Princeton  University  in  July,  1990.  He  arrived  at 

Princeton  after  9  years  as  Director  of  the  Aeronautical  Research  Laboratory, 
Department  of  Defense,  Melbourne,  Australia.  Prior  to  that  he  was  a  Professor  of 

Aeronautics  at  the  Graduate  Aeronautical  Laboratories  at  California  Institute  of 

Technology  and  a  senior  reader  and  lecturer  in  Mechanical  Engineering  at  the 
University  of  Adelaide,  Australia.  Prof.  Brown  has  published  widely  in  the  field  of 
turbulence  and  more  recently  on  vortex  breakdown.  His  work  with  Prof.  Roshko 
which  reported  the  discovery  of  coherent  structure  in  turbulent  mixing  layers  has 
been  recognized  by  the  Journal  of  Fluid  Mechanics  as  a  “classic  paper.”  His  work 
with  Wallace  (at  Adelaide  University)  and  subsequently  with  Mungal  (at  CalTech) 
developed  a  new  way  of  conducting  combustion  experiments  at  high  Reynolds 
numbers  in  highly  reactive  gases.  His  collaboration  with  Liepmann  and  Nosenchuck 
on  the  control  of  boundary  layer  transition  has  been  recognized  as  an  important 
result  in  the  development  of  turbulence  control.  Recent  work  with  Lopez  (at  the 
Aeronautical  Research  Laboratories)  on  the  classic  problem  of  axisymmetric  vortex 
breakdown  has  led  to  a  new  necessary  criterion  for  vortex  breakdown  based  on  the 
development  of  negative  azimuthal  vorlicity.  Prof.  Brown  is  continuing  his  work  on 

vortex  dynamics,  turbulence  and  compressible  gas  dynamics  in  collaboration  with 
other  members  of  the  department. 


Sample  Publications: 

1.  G.L.  Brown  and  A.  Roshko,  “On  Density  Effects  and  Large  Structure  in  Turbulent 
Mixing  Layers,”  J.  Fluid  Mcch.  64,  pp.  775-816  (1974). 

2.  P.E.  Dimotakis  and  G.L.  Brown,  “The  Mixing  Layer  at  High  Reynolds  Number: 

Large  Structure  Dynamics  and  Entrainment,”  J.  Fluid  Mech.  78,  pp.  535-560 
(1976). 

3.  G.L.  Brown  and  Andrew  S.W.  Thomas,  “Large  Structure  in  a  Turbulent  Boundary 
Layer,”  Phys.  Fluids,  20,  No.  10,  Pt.  II,  pp.  243-252  (1977). 

4.  H.W.  Liepmann,  G.L.  Brown,  and  D.M.  Nosenchuck,  “Control  of  Laminar- 
Instability  Waves  Using  a  New  Technique,”  J.  Fluid  Mech.  118,  pp.  187-200  (1982). 

5.  G.L.  Brown  and  J.M.  Lopez,  “Axisymmetric  vortex  breakdown,  Part  2.  Physical 
mechanisms,”  J.  Fluid  Mech.  (1990),  vol.  221,  pp.  553-576. 
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Anatol  Roshko 

Aeronautics  Department  105-50 
California  Institute  of  Technology 
Pasadena,  CA  91125 
(818)  356-4484 


Date  and  Place  of  Birth: 


[Pll  Redacted] 


B.Sc.  Engineering  Physics,  University  of  Alberta  1945 

M.S.  Aeronautics,  California  Institute  of  Technology  1947 

Ph.D.  Aeronautics,  California  Institute  of  Technology  1952 

Academic  Experience 

Instructor  (Mathematics)  University  of  Alberta  1945-46 

Lecturer  (Engineering)  University  of  Alberta  1949-50 

Research  Fellow  (Aeronautics)  Caltech  1952-55 

Assistant  Professor  (Aeronautics)  Caltech  1955-58 

Associate  Professor  (Aeronautics)  Caltech  1958-62 

Professor  (Aeronautics)  Caltech  1962- 

Theodore  von  Karman  Professor  of  Aeronautics  1985- 

Acting  Director,  Graduate  Aeronautical  Laboratories,  Caltech  1985-87 


Consultant 

McDonnell-Douglas  Corporation 
Rockwell  International 
General  Motors  Corporation 

Professional  Societies 

U.S.  National  Academy  of  Engineering 
American  Academy  of  Arts  and  Sciences  (Fellow) 

Canadian  Aeronautics  and  Space  Institute  (Fellow) 

International  Academy  of  Astronautics  (Corresponding  Member) 
American  Physical  Society  (Fellow) 

American  Institute  of  Aeronautics  and  Astronautics  (Fellow) 
Artie  Institute  of  North  America 
American  Association  of  University  Professors 
Sigma  Xi 

Wind  Engineering  Research,  Inc.  (Founding  Director) 

Fields  of  Interest 

(1)  Turbulent  Shear  Flow 

(2)  Separated  Flow 

(3)  Transonic  and  Supersonic  Aerodynamics 

(4)  Industrial  Aerodynamics 

(5)  Effect  of  Wind  and  Ocean  Currents  on  Structures 


Education 
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Page  2 


Committees,  etc. 


Associate  Editor,  Journal  of  Fluids  and  Structures  1986- 

Technical  Committee  on  Turbulence  Engineering  1987- 

Mechanics  Division  of  ASCE 

Advisory  Board  Stanford/Ames  Center  for  Turbulence  Research 
National  Research  Council  Aeronautics  and  Space  Engineering  Board  1988- 

Associate  Editor,  Physics  of  Fluids  1971-73 

Scientific  Liaison  Officer,  Office  of  Naval  Research,  London  1961-62 

National  Science  Foundation  Exchange  Scientist  to  India  1969 

U.S.  Wind  Engineering  Research  Council  Board  of  Directors  1978-85 

Member  National  Committee  for  Theoretical  and  Applied  Mechanics  1978-80 

Fluid  Dynamics  Panel  of  AG ARD  1984-89 

Honors 

Dryden  Research  Lecture,  AIAA  1976 

Turnbull  Lecture,  Canadian  Aeronautics  and  Space  Institute  1980 

Fluid  Dynamics  Prize,  American  Physical  Society  1987 

Division  of  Fluid  Dynamics 

Award  for  Professional  Achievement,  University  of  Alberta  1988 

Opening  Lecturer,  XVIIIth  International  Congress  of  1992 

Theoretical  and  Applied  Mechanics,  Haifa 
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External 

University  of  Houston  (Attn:  A.  Hussain) 

Princeton  University  (Attn:  S.  Orszag,  D.  Nosenchuck,  G.  Brown) 
California  Institute  of  Technology  (Attn:  A.  Roshko) 

University  of  Rhode  Island  (Attn:  F.  White) 

Massachusetts  Institute  of  Technology  (Attn:  K.  Breuer,  M.  Landahl) 
Yale  University  (Attn:  K.  Sreenivasan) 

University  of  Notre  Dame  (Attn:  M.  Gad-el  Hak) 

Brown  University  (Attn:  J.  Liu,  L.  Sirovich) 

University  of  Texas,  Austin  (Attn:  R.  Panton) 

Woods  Hole  Oceanographic  Institution 
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8233  (Myers,  Rice,  Varghese,  Krol,  Huyer,  Grant) 
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